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The density fluctuations that we observe in the universe today are thought to originate from 
quantum fluctuations produced during a phase of the early universe called inflation. By evolving a 
wavefunction describing two coupled Fourier modes of a scalar field forward through an inflationary 
epoch, we demonstrate that non-linear effects can result in a generation of entanglement entropy 
between modes with different momenta in a scalar field during the inflationary period when just 
one of the modes is observed. Through this mechanism, the field would experience decoherence and 
appear more like a classical distribution today; however the mechanism is not sufficiently efficient to 
explain classicality. We find that the amount of entanglement entropy generated scales roughly as 
a power law S cx A 1,75 , where A is the coupling coefficient of the non-linear potential term. We also 
investigate how the entanglement entropy scales with the duration of inflation and compare various 
entanglement measures from the literature with the von Neumann entropy. This demonstration 
explicitly follows particle creation and interactions between modes; consequently, the mechanism 
contributing to the generation of the von Neumann entropy can be easily seen. 

PACS numbers: 98.80.Cq, 98.80.Qc 
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I. INTRODUCTION 

Most modern cosmological models include a period in 
the universe's history called inflation during which the 
scale parameter increased exponentially with the proper 
time of a comoving observer. This period was originally 
introduced to address the horizon and flatness problems 
of cosmology More recently, however, research on 
inflation has been toward understanding structure for- 
mation @, [3, 4]. The distribution of galaxies and clusters 
that we observe in the universe today are thought to have 
originated from fluctuations of a quantized field created 
during inflation 0, 0]. A thorough review of structure 
formation and inflationary cosmology can be found in 
Liddle and Lyth Q . 

Despite their quantum mechanical origins, the late- 
time evolution of these fluctuations is treated in a clas- 
sical framework. It is therefore important to understand 
the quantum-to-classical transition made by these fluc- 
tuations (for a recent review, see Q). The classicality of 
a quantum system is often discussed in the context of de- 
coherence. That is, as a quantum system interacts with 
unobserved environmental influence, that system loses 
quantum coherence and begins to behave as a classical 
statistical distribution. 

The quantized field may of course be the inflaton it- 
self, which drives the inflation of the universe, or it could 
be another quantized field that produces density fluctua- 
tions as in curvaton models or the gravitational field. It is 
possible, in principle, that non-classical correlations from 
an inflationary period in our universe's history may one 
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day be observed. But this depends on the decoherence 
that the scalar or tensor field has experienced since the 
beginning of inflation. Several authors have investigated 
decoherence of the density fluctuations by calculating the 
entropy of cosmological perturbations created during in- 
flation 0, m m HEUE p. 

It has been suggested that decoherence is unlikely 
to occur during inflation because the Bunch-Davies state 
occupied by the scalar field during inflation is similar to 
the Minkowski vacuum. Because the ordinary Minkowski 
vacuum does not decohere, we would not expect to see 
any decoherence from a scalar field during inflation. In 
the particle-based picture adopted for the present anal- 
ysis, it becomes clear that the scalar field does undergo 
decoherence when the potential is non-linear. 

Since decoherence is a necessary condition for the 
emergence of classicality in a quantum system (T(| , non- 
linearities in the scalar field help to explain the classical 
matter distribution that we observe today. This sim- 
ple model demonstrates that this entropy generation can 
occur during inflation itself and does not depend on the 
reheating process at the end of inflation; therefore, the re- 
sults are perhaps most interesting for cosmological scalar 
fields that do not participate in reheating. For such fields, 
the non-linear interactions do not generate a sufficient 
amount of decoherence to result in classicality for the 
fields. 

Here, we examine the case where certain modes of a 
field play the role of the environmental influence and 
cause decoherence when a non-linearity in the potential 
allows the modes to interact [TH, [l?], EH ■ We discuss a 
simulation that was performed to compute the entangle- 
ment entropy between such modes in a very transparent 
model that follows particle creation and the interaction 
between modes during the inflationary period. The en- 



tropy is computed as inflation progresses to demonstrate 
the decoherence of a scalar field. 

Computing the entanglement entropy of a large quan- 
tum system is a computationally difficult task since it 
involves diagonalizing the density matrix. To evaluate 
several possible expediencies, we have compared our re- 
sults to other measures of entanglement and correlations 
between modes. We have found that the other measures 
considered share a similar qualitative behaviour with the 
entanglement entropy and can be much easier to com- 
pute. Therefore, for some applications, these measures 
may be useful as stand-in quantities in simulations where 
the entanglement entropy is too costly to compute. We 
verify several efficient methods to characterize the en- 
tropy. 



II. COSMOLOGICAL SCALAR-FIELD 
EVOLUTION 

We would like to investigate the evolution of a scalar 
field in an isotropic, homogeneous, flat spacetimc. The 
analysis for this situation is covered extensively in part 
I, chapter 6 of Mukhanov et al. . The relevant metric 
for this evolution is 



field. Even if the field itself is free, its self-gravity will 
introduce an interaction potential of the form [TBI, [l8| 



V = \M Pl (j> 3 . 



(3) 



Although the 3 potential is generally unstable, one 
should interpret this as an effective potential to account 
for the gravitational self-interaction, so the instability is 
not surprising because the gravitational self-interaction 
is generally unstable. 



A. Mode coupling during inflation 

For this analysis, we choose to use a simple model in 
which the universe contains only particles with four pos- 
sible momenta: ±k and ±2k. Given this requirement, 
we construct a Hamiltonian which incorporates a cou- 
pling term between these two Fourier modes so that we 
can observe the effect this non-linearity has on the en- 
tanglement between modes during inflation. 

The creation and annihilation operators satisfy the fol- 
lowing commutation relations 



[ok,4,] = <$( 3 >(k-k') 



(4) 



ds 2 



a 2 (r)(dr 2 -dx 2 ). 



(1) 



where r is the conformal time, which is related to the 
comoving time by dt = a(r)dr, and x is a comoving dis- 
placement. For simplicity we will take a(r) = — (Ht)~ 1 
(pure deSitter expansion) during inflation. 

The evolution of a scalar field <j> is governed by its La- 
grangian C. The lowest-order Lorentz-invariant expres- 
sion containing up to first derivatives is 



(2) 



For simplicity we will neglect the mass of the scalar field 
during inflation (m <C H). We include a non-linearity 
in the potential that couples the Fourier modes of the 



I a L a k'] = [ a k,a k '] = 0. 



(5) 



Including our potential term ([3]), the action for the 
field is 
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f 



d X\J — g 



d^d^cj) + XMp\(j) 3 
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Following the steps outlined by ref. [20(, we arrive at 
the following expression for the action. 



S = ^ / d 4 xa 2 



dcf> 



dr 



- [V0 2 ] + a 2 XM P ^ 



(7) 



J 



If we make the substitution u = a<j> = (2 J\a/2 J d 3 fcufc(r)e lk r , the action becomes 



= ^ / ird 3 k 



du k 



(k 2 + m 2 cS )\u k \ 2 - [ d 3 k'd 3 k"u^u^\\i + ^ + k") 



(8) 



where the effective mass is m 2 ff 



Q = 



and w = p/ p is the equation of state parameter. 
The Hamiltonian is, then, 



[(l-3w)} 



(9) 
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In general, we have itk = g(k,r)a] fL + 5*(fc,r)a_ k . Putting this into to the Hamiltonian, (jlOp . and neglecting terms 
that do not conserve energy in flat spacctime gives 



H = \ 



j d 3 k 
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dg* 2 (k,r) 
AM P1 



+ \9(k, r)\ 2 {k 2 + m 2 cff )j (4a k + a k a k ) 
+ g* 2 (k, r)(k 2 + m c 2 ff )) al k a k + + g 2 (k, r)(k 2 + m c 2 ff ) 
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The mode function is normally chosen to be 



g(k, r)g(k, r)g* (2k, ^a^a^a^ + g*(k, r)g* (k, r)g(2k, r)a 2 ka k a f 



(11) 



1 , , ,e~ lk7 



g(k,T) = —=(i-kr) (12) 

V2fc 3 r 

as this choice satisfies the equation of motion for the free field during a deSitter phase and because it simplifies the 
Hamiltonian to one that commutes with the number operator since, when Q = 1, 



dg\k,r) 
dr 



+ g'(k,T)(k'+m z eS ) = 0. (13) 



However, this choice is not practical for our calculation because the scalar field is not free; therefore, this choice does 
not satisfy the field equation of motion, and in fact it complicates the Hamiltonian because, for example 

g(k,T)g(k,r)g*(2k,T) = ^± J -(2(krf - 3t(kr) 2 - i) (14) 

does not have a simple dependence on r and the simplifications provided by (fl2"|) are lost. 

We would like to know the amount of entropy at the end of inflation during radiation domination. The usual way 
to proceed is to select the mode function (fT2|) and use this to determine the equation of motion for the scalar field 
during inflation. We would then determine the Bogoliubov coefficients at the transition from inflation to radiation 
domination. After performing the transformation, we would compute the amount of entropy from the transformed 
density matrix. 

However, we can simplify the problem by instead choosing a mode function that describes the system during 
radiation domination and use this mode function to compute the entire evolution. The choice of function g{k,r) is 
flexible due to the vacuum ambiguity and is related to choosing the set of states that the creation and annihilation 
operators act upon. Any choice will provide us with a complete basis with which we can describe any state of the 
field. The arrbitrariness of the mode function is also discussed in plj . 

For us, it is most prudent to choose the simple function 

g(k, T ) = -j=e- ikT (15) 

which defines the vacuum both during radiation domination and for scales much smaller than the horizon even during 
the de Sitter phase. Thus, we can make a very natural connection between our initial state and our final state. The 
choice is as arbitrary as choosing to perform a calculation in classical mechanics in a rotating frame rather than an 
inertial frame. 

Correctly interpreting the wavefunction where (|15|) is inappropriate (i.e. after horizon exit during a de Sitter phase) 
would require a Bogoliubov transformation, but for our purposes we do not require this. We are only interested in 
calculating the entropy after the transition to radiation domination where our choice of mode function corresponds to 
the usual creation and annihilation operators for this background. Therefore, we avoid transformations entirely since 
we already have the required description of our wavefunction. 

With the choice (|15p . the Hamiltonian is not constant in time even without the non- linear couplings. In particular 
the mass depends on time; this choice is similar in spirit to the calculations of Guth and Pi (Q). Heyl [2(| has shown 
that for a free scalar field that this choice gives the same results as the standard function g(k,r) and we refer the 
reader to that article for a more thorough discussion of the technique. 

Choosing to use (fT5"|) , we have 

Uk = -l=(e- ifeT a k + e ifcT al k ). (16) 



The nonlinear terms in the Hamiltonian provide a coupling mechanism between the modes of interest. To perform the 
integral over d 3 k' in (fT0|) . we neglect the effect of the coupling on the modes that are not considered in our simulation 
and treat the functions «k as constant on a spherical shell surrounding the momenta, k', that we are interested in. 
For Mk' = const, on spherical shells of constant volume around k and 2k, the integral becomes 



2k 



(17) 



where V — |7r ^Y+^f) ~ 23 is a (somewhat arbitrary) geometrical constant. 
Making this substitution, we arrive at the final form of the Hamiltonian. 



H 



d 3 k 



Q 



r 2 k 



Q 



-2ikr 



t t 2^kT^ Wk 3 / 2 M Pl . t t + . 

,T „t„2«a.t) + - — (4 k a k ak + a2k<a k ) 



+ a -k< e 



4V27ra 



(18) 



This Hamiltonian is similar to that used by ref. [201 ] . generalized to allow for the interactions between Fourier modes. 

Here, the two terms multiplied by the factor A are responsible for the annihilation of two particles from the k 
mode into a single particle from the 2k mode and the decay of an 2k mode particle into two k mode particles, 
respectively. As the two modes of the field exchange particles with each other, we expect that entanglement entropy 
will be generated in either of the modes observed individually. 

We wish to use this Hamiltonian to evolve Fock space wavefunctions representing the number of particles in each 
of four modes: Those with m + particles with momentum 2k', mT particles with momentum —2k', n + particles with 
momentum k', and n~ particles with momentum — k'. 
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Y, S m+ , m - in+in -(r)|m+,2k , )|m-,-2k')|n + ,k / )|n-,-k'). 

m+,m~,n+,n — — 

Whenever possible, we will use simplified notation such as 

oo 

\i>)= X B « 



(19) 
(20) 



: (r)|m ± )|n d 



(21) 



n^.m^— 



In order to evolve the wavefunction forward in time, we replace t with a new variable, x = —1/ (fcr). The equation 
of motion is then found from i-4-\ip) = H\i(}), left multiplied by (n^, ±k| {m il 1 ±2k|. The following identities are 
needed to evaluate H\tp): 



a k a kl TO± / 
(a k a k + a k a k 

d-kOklm, 



aL k a k |m ± ; 



n ) 

)IV>) 

n ± ) 



[m+S {3) (k - 2k') + n+6 {3) (k - k') + mS® (k + 2k') 
(2a k a k + Z)|V) 

Vm+m-(J (3) (k - 2k') + 5 {3) (k + 2k'))|m ± - 1)1^) 



-^ 3 )(k + k') 



|m ± )|n ± ) 



+Vn+n-((5 (3) (k - k') + 5^{k + k'))|m ± )|n d 



(22) 
(23) 

(24) 



: ) = y/(n+ + l)(n- + l)(o (3) (k - k') + S^(k + k'))^) \n± + 1) 



°2ka k a k \ m± ) \ n± ) 


+ \/(m+ + 


l)(m~ -+ 


- l)05 (3) (k - 2k') + 5 {3 \k + 2k'))|m ± + 1)1^) 


(25) 


= \J m+(n+ - 


f l)(n+ - 


h2)^(2k- 2k')|m+ - l)|m-)|n + + 2)|n - } 




a 2 k a k fl k k re± ) \ n± ) 


+^/m-(n~ 


- + l)(n- 


+ 2)<5 (3) (2k + 2k') |m+) \mT - 1) \n~ + 2) 


(26) 


= y/n+(n+ - 


l)(m+ - 


h 1)<5 (3) (k - k')|m+ + 1) \m-) \n+ - 2) \n~) 






+y/n~(n~ 


- l)(m- 


+ l),5( 3 )(k + k')|m+)|m- + l)\n+)\n- - 2) 


(27) 



where Z = [at, ajj = S^(k — k) is an infinite constant. 



After some algebra, we find the time evolution of the states is given by 
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l -^A m ± n ± (x) 



-2i*(/x 
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where the matrices A and B are related by a phase transformation 



A, 



,{x) 



-i(m + -\-7n~ -\-n + +n~ -\-Z)(j(x)-l) /x r> 



(■>■) 



(28) 



(29) 



with 7(x) = 2 + Qx 2 and Z is an infinite constant (related to the renormalization of the vacuum energy). The 
dimensionless constant a has the value . To arrive at equation (|28[) , we have ignored terms that involve modes 

dbik and ±4k since we are not concerned with how these modes evolve for our present purposes. 

I 



We begin the simulation for small values of x, well 
before the modes cross outside the Hubble length. At 
such a time, there has been a negligible amount particle 
production, so our initial wavefunction is simply the Fock 
vacuum, \tp)i = |m ± = 0)|n ± = 0). In the limit of - <C H 
ors < 1, this initial condition corresponds to the Bunch- 
Davies vacuum. During vacuum-energy-domination, the 
equation of state parameter, w, equals —1. Therefore, 
neglecting the mass of the scalar field, the value of Q is 
unity. 



B. Entanglement measures 

Discussions of dccohcrence rely on the notion of an en- 
vironment: a collection of degrees of freedom that inter- 
acts and becomes entangled with the system of interest. 
Our model is naturally separated into modes with differ- 
ent magnitudes of momentum. Noting that the entan- 
glement entropy does not depend on our choice of which 
set of modes is the environment and which is the system, 
we identify the modes with momentum ±2k with the 
environmental degrees of freedom and the modes with 
momentum ±k to be the system. 

This choice represents an entanglement due to coarse 
graining the internal degrees of freedom of the scalar field 
based on scale. One can think of the coarse graining as 
either being due to practical limitations in the observa- 
tions that can be made or as physical limitations such as 
a mode being entangled with a mode with a wavelength 
greater than the horizon size. The latter case is discussed 
in El. 



We measure the entanglement between modes using 
two different entanglement measures. The first of these is 
the entanglement or von Neumann entropy. The other is 
the linear entropy While the former is more common, the 
latter is easier to compute and scales monotonically with 
the entanglement entropy. Figure [1] shows a comparison 
between these two measures for a = 0.2. 

The density matrix of the above described system is 



p = WW 



B 1 



(30) 

|m ± )|n ± )<n ,± Km ,± | (31) 



i it — . n 31 ,m' 



and we assume that the modes with momentum 2k are 
inaccessible to measurement. This gives rise to a reduced 
density matrix obtained from tracing over the unobserved 
degrees of freedom. 



p N = Tr M p = 



i n± =0 
/ oc 



|V>><# 



(32) 



= E E E s m± , n± < ±>n , ± 1^)^1(33) 

n'±=0n±=0 \m±=0 / 

The von Neumann entropy is then a measure of the 
entanglement between the N system and the unobserved 
M system. 



N 



s 



-Tr(pArlnpjv) 



E^ ln ^ 



(34) 



i=l 



where the p^s are the eigenvalues of the reduced density 
matrix, p^ ■ A system with a finite Hilbcrt space spanned 



by N basis states will have a maximum entropy S n 
InN. 




and n' labels the Fock states. 
n + + n~ , each n' state is m - 
partition function can be written 



Since the energy is m — 
- 1 times degenerate, the 



/3m 



n' — l m—0 

Using the relation 



(e-/ 3 - l) 2 



(38) 



= J2 n P"& = E < n + l)e _/3 "(l " e- p ) 2 (39) 



n'=0 



we can eliminate /? for (n) using 



(40) 
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FIG. 1: The evolution of entanglement entropy, S, and lin- 
ear entropy, Sl = |1 — Tr(p 2 )|, as the mode stretches past 
the horizon for a — 0.2. This demonstrates that the non- 
linearities in the inflaton potential are capable of producing 
entanglement entropy between the coupled modes. Also note 
that S scales monotonically with Sl ■ 

The linear entropy, Sl = 1 — Tr(p 2 ), is often used as 
a stand-in for the entanglement entropy since it can be 
computed more easily and in our case contains the same 
qualitative information, 



i=i j=i 




(35) 



A system with a finite Hilbert space spanned by N basis 
states will have a maximum linear entropy Sl max = (N— 
1)/N. 

From figure [TJ we can see that this quantity is nearly 
proportional to the entropy. We will present the results 
both in terms of entanglement entropy and Sl ■ 



C. Thermal Entropy and Classicality 

The amount of entropy generated can be compared to 
the entropy of a thermal system that contains the same 
average number of particles. For a thermal system, the 
entropy is 

oo 

Sth = - e P n > th ln P n < th ( 36 ) 

71=1 

where the thermal density matrix is given by 



Eoc 
n> = l e 



-0E„ 



(37) 



where (n) is the average number of particles in the re- 
duced system. Finally, we can write the thermal entropy 
as 



Sth«n)) 



E 

771=0 



(m+1) 



4(n) T 



In 



(2+(n))™+ 2 V.(2 + (n)) 



4(n) 



\ 771+2 



(41) 

This quantity allows us to compare the entropy generated 
due to the coupling with the total energy of a thermal 
system at the same temperature. For example, if the 
information content of a system is defined as / = Sth — S 
then the relative information lost from the system due to 
the non-linear coupling term is 



-Aost — 1 — 



S 



S 



(42) 



Figure [2] shows that the rate of information loss due to 
the coupling is roughly the same as the rate of particle 
production. 

Campo and Parcntani argue that for Gaussian states 
at the threshold of separability and for (n) ^> 1, the 
entanglement entropy between modes will be one half 
the entropy of the thermal state [l3[. Since our states 
are not Gaussian, there is no known general separabil- 
ity condition. However, from the lack of growth in the 
information loss function shown in Figure [21 we can see 
that the Gaussian separability condition is unlikely to 
occur as (n) grows much larger than 1 at times greater 
than can be shown on the figure. Therefore, these types 
of non-linear interactions alone are likely insufficient to 
cause the system to appear classical. 

Another measure of separability used by Campo and 
Parcntani is the parameter S defined by the equation 



|c| 2 = n(n + 1 - S) 



(43) 



where n = ^(pa^a^) = and c = Tr (pa^a-^). The 
parameter, 5, is a measure of the correlations between 
the k and — k modes. For Gaussian density matrices, it 
can be shown that separability occurs when (5 = 1. The 




FIG. 2: The fraction of information lost due to tracing out the 
unobserved degrees of freedom, defined by equation (|42[) . and 
the separability parameter 5 (Eq. 143 1 as the observed mode 
stretches past the horizon for a = 0.2. By the end of the 
simulation, Ji os t appears to have leveled off to a constant few 
percent. For a separable (classical) system, we would expect 
host to grow to at least 0.5 to the right of the graph. If the 
system were to become classical, we would expect 8 to grow 
to 1. However, it too levels off to less than a percent. 



value of 5 measured for our model is shown alongside the 
information loss function in Figure O In both cases, the 
measures flatten out after the modes leave the horizon 
and fail to grow as one would need for non-linearities to 
explain the classicality of the quantum state. We can 
generalize the definition of c to measure the correlation 
between modes of different magnitudes of momenta in 
our system 



d = Tr(pa_2kakak) 



(44) 



Although the interpretation of this quantity or S is not as 
clear cut as for Gaussian density matrices, we find that 
both are useful and convenient tracers of the entangle- 
ment entropy. 



D. Estimating the sizes of A and £fl na i 

In order to match our above analysis with reality, we 
would like to make order of magnitude estimates for the 
parameters a in equation (|18|) and the final value of the 
x at the end of inflation, 2fi na i. 

For fluctuations in a scalar field other than the infla- 
ton, the value of A is essentially arbitrary; however, the 
gravitational self-interaction of the field provides a strict 
lower bound. Burgess, et al. [HI give an estimate of this 



self-interaction, 
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10 14 GeV 



(45) 
(46) 



where M 4 is the vacuum energy associated with the 

scalar field, and e = C^rJ i s a slow-roll parameter 
which may be larger than 1 if the scalar field is not the 
inflaton. We have included a possible matter-dominated 
period following the end of inflaton from scale factor oei 
to <zrh before reheating and taken a to be the value of 
scale factor at the end of inflation. 

The parameter a was introduced in equation l|18p to 
replace 



XV H 



a 



8y/2nkAI Pl 



(47) 



So, if we take, for example, a mode of size u> = ck ~ 
0.1Hz = 5 x 10~ 45 Mpi today, we arrive at an estimate 
for a due to gravitational self-interactions. 



2 x 10" 



0.01 



-3/2 



M 



10 14 GeV 



0.1 Hz 



-1/2 



(48) 

If the scalar field in question is the inflaton field, the grav- 
itational self-interaction will dominate over self-coupling 
interactions. 

The analysis here has assumed that reheating is quick 
and efficient 0, [23[, but in principle the end of infla- 
ton may be followed by a period of matter domination 
from scale factor <zei to orh before reheating. With this 
generalization, the comoving Hubble rate at the end of 
inflation is 
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30' 



7T a E i 
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Qgi 
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(49) 
(50) 



where M is the vacuum energy associated with the in- 
flaton field, (« AAfpj/4) and g r is the number of rela- 
tivistic degrees of freedom at the end of reheating where 
the photon counts as two. The value of Xfi na i (at the end 
of inflation) for the comoving scale a^H is simply unity 
and for other scales we have 



*£final 



a E iH 



6.3 x 10'^ 



7 l/4_ 



M 0.1 Hz 



10 14 GeV 



(51) 



Consequently although the correlations are present on 
all scales, they are most obvious on the comoving scale 
of the Hubble length at the end of inflation (i.e. really 
small scales). On these small scales the density fluctua- 
tions are well into the non-linear regime today but tensor 
fluctuations, gravitational waves (GW), would still be a 



loyal tracer of these correlations. Inflationary tensor per- 
turbations were first calculated in [24j]. 

The expression given in Eq. [51] is very uncertain. 
Typically today's Hubble scale is assumed to pass out 
through the Hubble length during inflation about 50—60 
e— foldings Q; Eq. ISTl gives 56 e— foldings before the end, 
so the centihertz scale would pass through the Hubble 
length 12—22 e— foldings before the end of inflaton. How- 
ever, the former number is highly uncertain. For exam- 
ple, if inflation occurs at a lower energy scale or if there is 
a epoch of late "thermal inflaton" [25|, [26|, [27J , the num- 
ber of e— foldings for today's Hubble scale could be as 
low as 25 0, yielding x nna i -C 1 for ck ~ 0.1 Hz. 

Because the simulation increases in complexity as par- 
ticles are produced (see figure|3]), we are confined to keep- 
ing xnnai ~ 0(1). So, even though a may be small in 
reality, there may be sufficient time during inflation for 
even a small non-linearity to produce a great deal of en- 
tanglement entropy because of very large values of Xfi na i ■ 



III. RESULTS 

We would like to investigate how the amount of en- 
tropy generated in a single mode scales with the coupling 
strength and the duration of inflation (i.e. a and xa na {). 
Figure Q] explicitly shows the creation of entanglement 
entropy for a = 1 as the universe undergoes its inflation- 
ary phase. The horizontal = H/p, is the physical 
size of a mode with respect to the horizon scale. The 
entanglement entropy increases less quickly than expo- 
nentially, which would be a straight line on the figure. 
Unfortunately, as was mentioned previously, the compu- 
tational size of the problem prevents us from simulating 
far past horizon crossing because the number of parti- 
cles becomes too large. Figure [3] shows how many Fock 
states are in the reduced system at each time step in the 
simulation. The number of states being integrated is this 
number to the 3/2 power, and the number of entries in 
the density matrix is the square of this number. 

The evolution of particles in the system is shown in 
figure [H Our results are consistent with those found in 
Heyl [20l | and show a nearly exponential evolution of the 
average particle number. Moreover, we can look at the 
evolution of each mode separately. For A = 0, each mode 
evolves according to the same equations of motion, and 
in this case, there is no difference between the rate that 
each of the modes evolves. However, the nature of the 
interaction between the modes is not symmetric because 
the decay of a single M mode particle results in 2 N 
mode particles and therefore the interaction results in an 
increased rate of production of N mode particles, relative 
to the M mode. Figure [5] shows how the entanglement 
entropy scales with average particle number when a = 
0.2. 

We performed the simulation for a variety of values for 
the coupling, a, spanning several orders of magnitude. 
Figure [HI shows entropy generation as a function of a for 
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FIG. 3: The number of Fock states associated with the re- 
duced system. 
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FIG. 4: Evolution of the average particle numbers for each 
mode for a — 1.0 

a variety of inflation durations x nna i- From this plot, we 
can see that S^nai scales roughly as a power law in a. 
Most of the a dependence can be removed by dividing 
Sfinai by a 175 . Doing this also helps to illustrate how 
Sfinai scales with x nna i- As expected, there is no entropy 
generated without the coupling terms (i.e. when a = 0). 
In this case, there is no communication between modes 
of the scalar field and they evolve independently. 

As was mentioned earlier, Sl is a useful stand-in for S 
that can be computed faster than S. Figure [8] echoes the 
previous results in terms of Sl instead of S. In this case, 
1 — Tr(p 2 ) scales more like a 2 instead of a 175 . How- 
ever, both Sl and S demonstrate the same qualitative 
behaviour. 

In addition, we have found other useful stand-ins for 
the entanglement entropy that are easier to compute and 
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FIG. 5: Entanglement entropy vs. average particle number 
for a = 0.2 
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FIG. 6: von Neumann entropy vs. a for various values of 
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scale similarly with a. The S parameter, defined in 
scales roughly like an a 2 power law much like Sl ■ Figure 
[5] shows the power law behaviour of this function. Addi- 
tionally, if we use a simple measure of correlation between 
entangled modes, \d\ 2 fEq.EPIj), we find that its scales like 
q.i.85 ^ see gg [yjj ^ anc j so can De a US eful stand-in for the 

von Neumann entropy, S oc a 175 . 

In the real universe, we are dealing with small values 
of A and very large values of x. However, the simulation 
outlined in this paper is limited because its computa- 
tional complexity increases dramatically as particles are 
produced, even for small values of the coupling, a. More- 
over, for small values of a, the production of entropy is 
too small to be meaningful. While the dependence of S 
on a nearly follows a power law, there is no simple rela- 
tion describing the dependence of S on a;g na i. The value 



FIG. 7: \d\ 2 = |Tr(pa2kdkCik)| 2 scales with a in much the 
same way as S, but is less costly to compute. 
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FIG. 8: Sl vs. a for various values of Xfinai 

Sl is approximately proportional to a 2 x 3 over a wide 
range of a and the modest range x probed by the simu- 
lations; therefore, very roughly, we can write the scaling 
law as S l oc a 2 a;| nal where 

s ^ 101 ^ /4 (lo4v) I5 (^)" 4 - < 52 > 

Of course, only values of Sl less than unity make sense, 
so a larger value from the fitting formula indicates that 
S l is very close to one. However, a value of Sl < 1 is 
obtained by lowering the mass scale of inflation below 

M < 8 x 10 12 (ttj^) ~ 4/1 ° GcV; (53) 

therefore, if the energy scale of inflation is low, the quan- 
tum states of fluctuations at w ~ 0.1 Hz will remain 
coherent despite the non-linear coupling. 
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FIG. 9: 5, defined in equation (|43[) scales with a in much the 
same way as Sl, but is less costly to compute. 

The simulation was checked for consistency in several 
ways. First, we traced the probability throughout the 
simulation measured both by the sum of squares of the 
matrix elements J3m± n ± =o ■^m ± »± an d the trace of the 
density operator. Both of these quantities were conserved 
to a few parts in 10 -7 . Moreover, we estimated the level 
of numerical error by rerunning the simulation with a va- 
riety of phase rotations multiplying the initial wavefunc- 
tion. The standard deviation of the results from these 
numerical changes in the initial conditions give us an idea 
of the level of numerical error in the simulation, which 
were typically at the level of one part per thousand. 

IV. CONCLUSIONS 

In this paper we have developed a model in which two 
modes of a scalar field evolve during inflation and we 
have computed the entanglement entropy between them. 
The entanglement entropy generated between observed 
and unobserved modes in the inflaton field give the ap- 
pearance that entropy is being produced, even though 
the scalar field remains in an overall pure state. The 
preceding results clearly show that non-linearities in the 
inflaton potential give rise to a generation of entangle- 
ment entropy between observed modes and unobserved 
modes in a scalar field during inflation. This entropy is 
an additional source to that caused by coupling to ex- 



ternal degrees of freedom [28| , entanglement between the 
inside and outside of the horizon [29( and that which is 
created during reheating after inflation has ended. 

We have attempted to extrapolate the results of our 
simulation to the real universe. The relevant parameters 
determining the amount of entropy generated via non- 
linearities are the strength of the coupling a g ~ 10 -3 and 
the scale of the fluctuation at the end of inflation given 
by the dimensionless parameter Zfinai ~ 10 7 . The entan- 
glement entropy was found to scale like a 1 ' 75 for a fixed 
affinal ■ The dependence of Sl on Xfinai for a given value of 
a is not as straightforward, but Sl °c x^ nai over a short 
range of Xfmai values. Based on these rough scaling pat- 
terns, we estimate that non-linearities due to gravity and 
inflaton self-coupling are insufficient to decohere modes 
that spend only a few Hubble times at super-horizon 
scales. In particular, if the energy scale of inflaton is 
less than 10 13 GeV, fluctuations at about 0.1 Hz may 
remain coherent. 

We found two measures of the decoherence related to 
the correlations between modes of different momenta pro- 
vide a faithful estimate of the entanglement entropy in 
our model — one of these measures is new to this work 
(d) and specifically probes the non-linear coupling be- 
tween modes. In particular these estimates are very in- 
expensive to calculate as compared to the von Neumann 
entropy and should prove useful for more detailed models 
of entropy generation. 

It is usually assumed that the main contribution to the 
entropy observed in the density perturbations is gener- 
ated during reheating, when the inflaton decays. How- 
ever, the analysis demonstrates that entropy can be gen- 
erated independently of reheating provided there is even 
a small non-linearity in the scalar potential; therefore, 
the results are applicable to scalar fields that do not par- 
ticipate in reheating. For example, the gravitational wave 
background can be treated as a pair of scalar fields, so 
even tensor fluctuations may contribute to the entropy 
and the classicality of the distribution of density per- 
turbations in this way and observations of the gravita- 
tional wave background at high frequency could reveal 
the quantum mechanical origin of density fluctuations. 
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